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摘 要 :研究 了 Johnson 提出 的 倾 转 旋翼 不 平衡 载荷 前 飞 动 力学 模型 ,将 其 奖 叶 分 析 方法 应 用 于 直 
升 机 旋 收 系统 模 态 分 析 。 在 刚性 条 件 假 设 下 推导 了 直升机 旋翼 弹性 阻尼 和 惯性 力 综合 作用 时 桨 叶 
的 挥 攻 和 摆 振 运动 方程 ,给 出 了 固定 和 旋转 坐标 系 下 对 应 的 运动 方程 。 通 过 引入 均匀 入 流 和 线性 
扭转 假设 ,获得 了 运动 方程 的 理论 解析 解 。 利 用 登 加 原理 ,得 到 了 桨 载 轴 心 运 动 方程 ;采用 New- 
mark 法 进行 振动 微分 方程 求解 ,最 终 得 到 了 直升机 旋翼 的 轴 心 运动 轨迹 。 以 某 型 直升机 旋翼 系统 
为 例 ,验证 了 本 研究 所 提出 旋翼 桨 时 模 态 分 析 方法 的 准确 性 ,给 出 了 兼顾 计算 精度 和 效率 的 最 佳 求 
解 步 长 选取 方法 ;预测 了 典型 飞行 状态 下 的 桨 载 轴 心 运动 轨迹 ,为 直升机 旋翼 系统 设计 提供 了 基础 
方法 和 技术 参考 。 
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Abstract :This paper studies the forward flight dynamic model of tilt rotor unbalanced load proposed by 
Johnson ,and applies its blade analysis method to the modal analysis of helicopter rotor system. Under the 
assumption of rigidity ,the flapping and shimmy motion equations of helicopter rotor blade under the com- 
bined action of elastic damping and inertial force are derived ,and the corresponding motion equations in 
fixed and rotating coordinate systems are given. By introducing the assumption of uniform inflow and linear 
torsion ,the theoretical analytical solution of the equation of motion is obtained. Then ,using the superposi- 


tion principle ,the axial motion equation of the hub is obtained; Newmark method is used to solve the vibra- 
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tion differential equation ,and finally the axis motion trajectory of helicopter rotor is obtained. Taking a hel- 


icopter rotor system as an example ,the accuracy of the rotor blade modal analysis method proposed in this 


paper is verified ,and the optimal solution step selection method considering both calculation accuracy and 


efficiency is given. The hub axis trajectory under typical flight state is predicted, which provides a basic 


method and technical reference for the design of helicopter rotor system. 
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由 于 直升机 旋翼 能 够 同时 提供 升力 、 操 纵 力 和 
推进 力 , 使 得 直升机 具有 垂直 起 降 、 机 动 性 强 、 适 用 
范围 广 等 性 能 和 优点 ,但 是 同时 也 存在 着 其 特有 的 
飞行 稳定 性 问题 "2 。 直 升 机 通过 旋翼 旋转 在 桨 盘 
上 下 平面 产生 压力 差 进而 产生 升力 ,在 前 飞 过 程 中 ， 
旋翼 的 气流 状态 为 平 飞 来 流 、 纵 向 压 差 气流 和 旋 恤 
涡流 的 释 加 ,在 奖 盘 平面 上 产生 复杂 的 流动 速度 分 
布 ,并 产生 不 均匀 且 周 期 性 变化 的 气动 力 。 在 直 升 
机 慌 行 状态 下 ,旋翼 桨 叶 将 围绕 一 个 固定 牌 直 轴 转 
琉 孩 得 称 为 桨 载 旋翼 转轴 ,其 旋转 中 心 称 为 旋 避 轴 
心 本 在 周期 性 变化 气动 载荷 的 影响 下 ,旋翼 桨 载 轴 
心 阳 偏 离 平 衡 点 2 。 若 桨 载 轴 心 的 运动 随时 间 收 
人 钱 : 其 振动 幅度 不 断 减 小 , 且 振 动 频 率 相对 稳定 , 则 
站 天 机 处 于 动 稳定 状态 ; 若 浆 坑 轴 心 偏离 平衡 点 且 
旋转 发 散 ,振动 幅度 不 断 加 大 ,将 导致 直升机 机 体 发 
符 强 烈 颤 振 , 进 而 影响 机 体 飞行 性 能 ,严重 时 会 酿 成 
机 毁 人 亡 的 重大 安全 事故 。 因 此 ,国内 外 学 者 对 旋 
翼 和 浆 叶 的 气 弹 耦 合 分 析 ,振动 和 颤 振 机 理 稳 定性 
和 和 动 响应 分 析 \ 减 振 优化 设计 、 振 动 试验 等 多 方面 开 
展 下 大 量 卓 有 成 效 的 研究 。Nixonc] 详细 阐述 了 倾 
转 诡 器 的 颤 振 机 理 , 建立 了 倾 转 旋 翼 复 合 材料 桨 叶 
的 仿 弹 耦合 模型 , 开展 了 桨 叶 气 弹 动 响应 和 旋翼 系 
统 稳定 性 分 析 。 王 浩 文 等 四 针对 桨 叶 结构 ,提出 了 
一 种 气动 /惯性 耦合 的 旋翼 系统 分 析 方法 ,给 出 了 桨 
叶 前 飞 状态 下 的 非 线 性 周期 时 变动 力学 方程 及 旋翼 
系统 动 响应 。 文 献 [5-6] 则 基于 隐 式 多 浆 叶 坐标 转 
换 方 法 ,建立 了 一 种 前 飞 状 态 直 升 机 旋翼 /机 体 耦 合 
动 稳定 性 分 析 模型 ;并 在 此 基础 上 ,研究 了 直升机 前 
飞 室 中 共振 时 各 自由 度 的 相互 作用 ,为 揭示 该 共振 
物理 本 质 英 定 了 坚实 基础 。 王 红 州 等 "1 针对 某 无 轴 
承 旋翼 ,利用 混合 优化 策略 开展 了 气 弹 动力 学 多 有 目 
标 优化 设计 , 减 振 效果 显著 。 文 献 [8-9] 则 分 别针 对 
螺旋 桨 奖 叶 和 直升机 模型 奖 叶 开展 了 模 态 和 颤 振 试 
验 分 析 , 给 出 了 一 些 规律 性 的 试验 结论 ,为 桨 叶 减 振 
设计 提供 了 有 力 技术 支撑 。 

在 系统 总 结 和 梳理 前 人 研究 成 果 的 基础 上 ,本 


研究 基于 Johnson "提出 的 四 自由 度 模 型 ,引入 均 
匀 入 流 和 线性 扭转 的 假设 ,获得 了 桨 叶 挥 舞 、 摆 振 运 
动 方程 ,利用 县 加 原理 ,得 到 了 桨 坑 轴 心 运动 方程 ， 
并 使 用 Newmark 法 进行 振动 微分 方程 求解 ,最 终 得 
到 了 直升机 旋 辟 的 轴 心 运动 轨迹 。 最 后 以 某 型 直 升 
机 旋 壳 系统 为 例 ,计算 了 各 个 模 态 的 响应 形式 和 固 
有 频率 ,并 对 照 试 验 结果 进行 了 计算 误差 评估 。 


1 直升机 不 平衡 旋 蝇 动力 学 建 模 


1.1 旋翼 气动 力 模型 


本 研究 采用 叶 素 理论 ” 对 旋翼 气动 力 进行 建 
模 , 将 其 分 解 为 1 个 叶 素 ( 二 维 避 型 ) ,依据 普 衣 特 升 
力 线 理论 , 叶 素 的 受 力 示 意 如 图 1 所 示 。 以 单独 旋 
辟 为 研究 对 象 时 ,通常 定义 2 个 坐标 系 :旋转 坐标 系 
和 固定 坐标 系 ,2 个 坐标 系 之 间 相 差 一 个 绕 转 轴 的 
欧 拉 角 ( 奖 叶 方位 角 ) 。 旋 转 坐 标 系 常用 于 描述 单 
片 桨 叶 的 力 关系 ,一 般 是 波动 的 ;而 固定 坐标 系 常 用 
于 描述 整个 浆 盘 的 升力 ,一般 是 稳定 的 。 


dr 1 


图 1 叶 素 受 力 示意 图 
Fig. 1 Force diagram of leaf element 


图 1 中 :dQ 和 d7 分 别 为 dR 在 旋 波 固定 坐标 系 


内 水 平和 垂直 方向 上 的 分 量 ,表示 单个 叶 素 产生 的 
升力 和 阻力 ;Bu 为 来 流 角 ;eo 为 奖 叶 安 装 角 ;ao 为 迎 
角 ;r 为 奖 叶 半 径 ;2 为 旋转 角速度 。 
旋翼 奖 叶 合 速度 为 
U= VUs+ Us, (1) 
式 中 UV 分 别 为 叶 素 弦 向 和 垂 向 速度 。 
定义 升力 L 和 阻力 D,9 为 操纵 角 ,p 为 浆 叶 安 


装 夹 角 , 则 气动 迎 角 a =0 - p, 总 气动 力 垂直 于 和 平 
行 于 桨 盘 平 面 的 分 量 分 别 为 已 和 成 。 升 力 和 阻力 
的 公式 为 


=3p UeC, (2) 
D=3p UeC, (3) 
其 中 :e 为 桨 叶 弦 长 ;Ci 为 缀 型 的 升力 系数 ;Cj 为 蛋 
型 的 阻力 系数 。 
气动 力 在 两 个 方向 上 分 力 的 大 小 分 别 表示 为 
F.=Leosp - Dsing (4) 
F,=Lsing + Deosp (5) 


浆 叶 上 作用 的 基 元 拉力 、 基 元 扭矩 及 基 元 功率 
可 表示 为 


dT =N,F.dr (6) 
dO =N,F,dr (7) 
dP =wdQ =w N,F,dr (8) 


式 中 ， ,为 桨 叶 个 数 ;r 为 桨 时 半径 ;io 为 桨 盘旋 转角 
带 度 ;P 为 功率 。 通 过 对 上 述 基 元 载荷 进行 积分 即 
可 得 到 旋 织 的 整体 作用 力 o 
1, > 旋翼 轴 心 运动 分 析 


局 在 刚性 条 条 件 下 , 桨 叶 的 挥舞 运动 方程 “可 写 为 
大 二 天 本 一 27 peBE + Ta, (a, Ra 一 


(过 = 2 Ge )cos pn， + Sg En = | qr = yM, 
全 (9) 
其 中 万 是 量 纲 归 一 化 旋翼 的 挥舞 惯性 矩 , 写 为 
全 I R 
1s = 了 1, = 上 mrdr (10) 


式 中 ， m 表示 桨 叶 微 段 质量 ;y 为 旋 避 洛克 数 ,数值 
等 于 气动 力 和 惯性 力 的 比值 ;vw 为 旋翼 的 挥舞 自然 
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E=é& + .cosny + ésinny) (13) 
其 中 :a 为 升力 线 系数 ;名 为 旋 芝 的 平均 摆 振 角 ,反映 
了 旋 避 上 的 合 扭矩 对 旋 必 的 作用 效果 ; 台 。 6 为 旋 可 
整体 的 偏 移 , 是 旋 枝 对 周期 性 扭矩 的 啊 应 。 一 般 只 


考虑 &1. 1.。 
是 ,在 固定 坐标 系 下 ,旋涡 挥舞 运动 方程 可 进 
一 步 写 为 


-2 有 BEo+S 2,=Yy Mr (14) 


I B(é +é,) - 


1 [6 + vgsBo] 
Ls [Bi +2B,, + (vs -1)Bi.] -2 


I (a, -2 60)=Yy Mn, (15) 
[By,-2B,+ (vs -1)B,] +2 BC 一 所) - 
Ls,( a, +2 oo ) =Yy 有 (16) 

同 理 , 在 固定 坐标 系 下 ,旋翼 摆 振 运动 方程 可 进 


一 步 写 为 

及 (E+zeto)-2 友 BBo -lia,=y Mo (17) 
le [E+2 6 + (ve -1)é.] -2 B(Bi, +B.) - 
la( Qa, -2 0) -Se y=Yy Ma (18) 
-2 + (m1)8,) +218p(B, -Bi) ~ 
及 (as +2 6) +S x, =Yy M,, (19) 
其 中 :as or Qs 和 wi yz 分 别 为 桨 载 轴 心 3 个 方向 
的 角 位 移 和 线 位 移 。 


2 直升机 旋翼 轴 心 运 


Li(é, 


云 动 轨迹 计算 


2.1 运动 方程 求解 过 程 
已 知 旋 弃 气动 力 的 情况 下 ,运动 方程 中 Bb, 等 稳 
态 数值 可 通过 迎 角 、 浆 盘 倾角 和 挥舞 角 的 几何 关系 


频率 ,对 于 匀 接 式 旋 翼 约 为 1/rev。 

相应 地 , 浆 叶 的 挥舞 运动 方程 在 旋转 坐标 系 下 
可 表达 为 

B=Bot+ + ,8, .Cosny + B,,sinny) (11) 

式 中 :Bo 为 稳 态 锥 度 角 ; B,, 和 BB, 为 桨 叶 叶 尖 平 面 偏 
离 挥 舞 平面 的 角度 ;y 为 桨 叶 方 位 角 。 由 于 高 阶 谐 
波 分 量 反 映 了 桨 尖 轨迹 平面 的 手 曲 模 态 , 故 描述 直 
升 机 前 飞 挥舞 运动 时 ,主要 用 一 阶 8;. 和 B,, 表 示 。 

同 理 , 在 刚性 和 和 柔性 条 件 下 , 旋 惨 摆 振 运动 方程 
分 别 可 表达 为 


,nF 
I (E+ vt) +2BB = y[ ndr (12) 


求 得 。 实 际 情形 中 , 随 着 频率 的 增 大 ,其 谐 波 越 来 越 
低 , 对 应 的 运动 幅度 也 越 来 越 小 ,只 有 前 若干 阶 运动 
有 实际 意义 ,而 后 续 的 表达 式 均 为 控制 桨 叶 的 弯曲 
模 态 数值 ”。 因 此 可 以 考虑 使 用 前 对 级 截断 的 方 
法 ,而 后 使 用 级 数 对 比 的 方法 求 得 8 和 B,.。 

桨 载 力矩 可 表示 为 


—2 Cy, 

Oa -De (20) 
2 Ci 4 Bi 

Oa 


其 中 :Cw 、Cw 为 对 应 的 力矩 系数 ;0 为 旋 惨 实 度 ;a 
为 升力 线 系数 ;y 为 旋 避 洛 元 数 。 通 过 式 (14) 可 将 
桨 叶 的 运动 模 态 和 受 力 状 态 联系 起 来 , 从 而 进 
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行 计算 。 

针对 贸 接 式 旋 渗 , 挥 舞 贸 的 刚性 转动 模 态 可 表 
达 为 


k(R-e,),R=e, 
9=10 Re 
式 中 e, 表 示 挥 舞 贸 偏 置 量 。 于 是 ,旋翼 的 拉力 系数 
表示 为 


(21) 


0,,， 1 
| 
Kp er 
$8 (2 
其 中 为 浆 叶 流入 比 。 另 一 方面 ,前 文 已 经 给 出 挥 
舞 运动 自然 基 阶 频率 的 表达 式 ,引入 质量 均匀 假设 ， 
其 结果 可 简化 为 
, 2 3e, 
> ”| 
好 可 将 挥舞 模 态 与 旋翼 所 受 的 拉力 相 联系 , 摆 振 运 
政 守 此 同 理 ,对 应 的 表达 式 分 别 为 


Ee Oa Oo 有 3 2 大 EE On 2 
| (La | 8 


(23) 


区 


1 e 
ea 0 Hp 
7 (A MK 01.) Br] 三 本 (24) 


CN 对 于 系数 周期 为 7 的 动力 系统 ,可 用 如 下 形式 
所 好 为 


X=A(i)x+B(i)u (25) 
-一 满足 条 件 
pe A(t+T) =A(t) (26) 
则 其 解 可 表达 为 


r= = fx 0+ [f(s) Bs)uds (27) 
“由 指数 函数 的 性 质 可 知 ,周期 系数 方程 必 存 在 
如 下 特 解 。 

fli,t0) =Pe" =0e"0- (28) 

这 里 矩阵 8 是 矩阵 P 了 分解 出 的 对 角 和 矩阵 。 

考虑 其 初始 条 件 , 令 X+、Xp 分 别 为 X'(0) =0， 

X(0)=1 和 X'(0) =1,X(0) =0 两 个 状况 下 的 特 
解 ,构造 特 征 和 矩阵 和 特征 方程 为 
寅 。 充 。 

| 5 

s — (Xa+Xp)s + (XX, XX) =0 (30) 

中 就 包含 了 对 强迫 输入 量 wu 的 振动 响应 , 即 可 获 

得 旋翼 在 气动 力作 用 下 的 各 模 态 响应 。 


2.2 模 态 数据 求解 方法 
求解 过 程 是 移 根 据 机 型 数据 ,计算 出 旋翼 在 不 
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同 状态 下 所 受到 的 总 体力 和 力 抢 ,计算 出 旋翼 桨 叶 
的 前 几 阶 固有 模 态 ;然后 采用 Newmark 法 于 ”进行 
迭代 求解 ,计算 得 到 旋翼 奖 叶 的 响应 ,将 其 三 加 得 到 
旋翼 轴 心 运动 轨迹 ; 使 用 相关 软件 输出 其 运动 图 
像 ,根据 图 像 的 具体 情况 来 直观 判断 运动 是 否 发 散 ， 
若 发 散 则 分 析 其 具体 发 散 模 态 ,和 若 未 发 散 , 则 说 明 旋 
辟 工 作 正常 。 


符 


用 Newmark 守 .让 笛 某 吉 种 心得 移 


给 出 柴 吉 轴 


Li 
村 前 禄 是 ， 得 到 运 引 方程 解 相 鱼 
到 2 ” 轴 心 运动 计算 流程 图 
Fig.2 Flow chart of axial motion calculation 
在 已 知 直升机 旋翼 质量 和 刚度 分 布 的 情况 下 ， 
可 采用 牛顿 法 有 限 元 法 等 求解 旋 咽 的 挥舞 、 摆 振 、 
扭转 等 固有 特性 ,为 下 一 步 轴 心 运动 轨迹 计算 提供 
输入 。 下 面 以 牛顿 法 为 例 , 介 绍 旋翼 固有 频率 计 
算 过 程 。 
旋 壳 挥舞 运动 的 控制 微分 方程 一 般 写 为 
(EJsy)"-(Ny')’+my =0 (31) 
其 中 :E J 为 桨 叶 刚 度 ;N 为 桨 叶 上 的 离心 力 ;m 为 
浆 叶 微 段 质量 ;7 为 染 叶 刚体 位 移 。 
此 方程 可 使 用 分 离 变量 法 求解 ,或 直接 借助 计 
算 机 工具 数值 求解 。 此 处 采用 分 离 变量 法 求解 ,将 
其 分 离 变 量 后 积分 两 次 ,可 得 
[Ey Cy) dr + [Nd — ,| mdr =0 
(32) 


分 村 发 散 悦 起 


若 第 i 阶 广义 刚度 和 广义 质量 写 为 
K = [EJ (%)dr + Cy)dr| mdn (33) 


R 
m= fidr (34) 
则 旋 中 的 第 i 阶 挥舞 固有 频率 为 
Kk. 
wp = (35) 


将 相关 数据 带 入 1. 2 节 挥 舞 摆 振 运动 控制 方 
程 ,可 以 获得 匀 接 和 无 匀 ( 柔 性 ) 式 旋 豆 的 边界 


条 件 为 
y(0)=0, JR) =0，Y 0) =0 (36) 
y(0)=0, 7y(0)=0, y'(R)=0 (37) 
容易 求 得 旋 波 挥舞 的 0 阶 固有 频率 为 


WM, 
om = /1 +L ) (38) 
已 


式 中 :万 为 桨 叶 的 铵 外 伸 量 ;7 为 挥舞 匀 轴 的 静 和 矩 ; 
11 为 挥舞 匀 轴 的 惯性 矩 。 

同 理 ,可 求 得 旋翼 摆 振 和 扭转 运动 控制 微分 方 
程 的 积分 表达 式 和 固有 频率 ,此 处 就 不 再 袭 述 。 


2.3” 轴 心 运动 计算 


本 研究 主要 利用 Newmark 法 “中 进行 响应 计 
算 。 在 具体 计算 时 ,关键 在 于 通过 目标 旋翼 的 运动 
阿 量 和 质量 .阻尼 刚度 矩阵 ,计算 出 稳 态 状态 下 桨 
根 办 矩 表达 式 , 令 w=0 可 得 到 桨 叶 在 0° 方位 角 时 
的 余力 初始 值 ,可 作为 最 初 的 力作 用 向 量 。 为 便于 
迭代 计算 ,初始 运动 位 置 可 设 为 位 于 预定 的 浆 载 坐 
标 原 点 ,在 时 间 单位 化 后 即 可 获得 在 某 个 特定 时 间 
或 特定 时 间 段 内 的 位 移 数 据 """。 旋 严 3 方向 质量 
列 度 和 矩阵 可 由 工程 梁 理 论 求 得 , 带 入 样板 直升机 相 
芜 质 量 和 惯性 矩 数 据 即 可 计算 ,相关 推导 如 下 。 
@@ 桨 叶 弯 和 矩 和 桨 叶 变 形 的 本 构 关系 式 为 
二 M.(7) -E19 (39) 
泛 对 应 的 弯曲 运动 微分 方程 为 


jd dl dr 加 
@ (| mpdp FE)- mds +t ms =F 


到 (40) 
吕 另 一 方面 ,将 前 文 得 到 的 各 模 态 数据 求 出 , 即 可 
求 得 桨 叶 整 体 的 运动 方程 ,进而 求 出 桨 载 轴 心 的 运 
动 轨迹 ,相应 的 坐标 变换 为 


X B Qt 0 x 0 
了 |=| < | 4 | X : (41) 
Z 0 OQ 丰 0 0 


将 上 述 各 式 联 立 ,以 矩阵 形式 输入 相关 机 型 数 
据 , 将 气动 力 数据 写成 向 量 形式 ,作为 前 文中 提 到 的 
u 输入 , 即 可 获得 其 计算 表达 式 。 


3 ”直升机 旋 辟 轴 心 运动 分 析 与 轨迹 预测 


3.1 旋 收 桨 叶 模 态 求解 分 析 
本 研究 中 选取 某 型 直升机 作为 对 象 进行 算 例 分 


析 , 相 关机 型 数据 如 下 表 。 
表 1 某 型 直升机 相关 参数 


Tab.1 Relevant parameters of a helicopter 


旋 旨 数量 3 片 
辟 型 NACA0012 
桨 叶 挥 舞 惯 性 矩 /(kg .ma ) 3 891. 198 
奖 叶 极 惯性 矩 /(kg'm) 15 727. 489 
洛克 数 8.1 
旋 经 实 度 0.085 
贸 外 伸 量 比 0.05 


机 身 三 方向 惯性 矩 /(kg . m2 ) 5423.272 6779.09 40674.542 


平 飞 叶 尖 速度 /(m : s-1) 198. 12 
浆 叶 弦 长 m 0.3 
浆 盘 半径 /m 5.345 
扭 度 /(°) —12.275 
扭转 方式 线性 扭转 


以 上 数据 代入 式 (39) ~ 式 (41) ,其 相关 模 态 固 
有 频率 数据 求解 如 下 。 
表 2 步 长 0.1 的 各 模 态 计算 结果 


Tab.2 Calculation results of each mode with steps 0.1 


模 态 名 称 实测 值 / 步 长 0.1 计算 误差 | 计算 时 间 / 
Hz 计算 值 /Hz % s 
一 阶 挥舞 15.75 16. 82 6.79 3.841 
一 阶 摆 振 19.51 21: 37 9.53 4.677 
一 阶 扭转 144.70 157.30 8.71 4.251 
二 阶 挥舞 41.93 45. 14 7.66 8.912 
二 阶 摆 振 54.68 58.99 7.82 10. 022 
二 阶 扭转 193. 86 202.03 7.88 9.736 
三 阶 挥舞 86.84 94.55 8.88 14.715 
三 阶 摆 振 97.74 103. 15 8.60 16.814 
三 阶 扭转 243. 12 250.76 5.20 16. 129 


依据 表 2 ~ 表 7 的 计算 数据 ,将 各 模 态 在 不 同步 
长 下 计算 所 得 相对 误差 及 所 需 计 算 时 间 分 别 作 图 
(图 3~ 图 4) ,为 了 方便 对 比 ,在 横 轴 绘制 时 未 使 用 
均匀 刻度 而 使 用 了 对 数 刻度 。 本 研究 6 次 计算 所 使 
用 的 步 长 分 别 为 :0. 1.0. 05 .0. 01 .0. 005 .0. 001、 
0. 000 5。 对 比 数据 发 现 ,计算 得 到 的 模 态 固有 频率 
比 实测 值 偏 大 ,主要 由 于 简化 条 件 引 起 的 系统 误差 。 
观察 图 3 和 图 4 可 知 : 随 着 计算 步 长 减 小 ,前 三 阶 模 
态 计 算 值 也 越 来 越 逼近 实测 值 , 但 计算 时 间 也 逐渐 
增加 ;在 计算 步 长 减 小 到 0.01 以 后 ,前 三 阶 模 态 计 
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算 误差 已 减 小 到 5% 以 下 , 且 计 算 误差 变化 很 小 。 


表 3 步 长 0.05 的 各 模 态 计算 结果 


Tab.3 Calculation results of each mode with steps 0.05 


并 起 名 称 “实测 值 ” 。 步 长 0.05 。 计算 误差 计算 时 间 / 
Hz 计算 值 /Hz % s 
15.75 16. 44 4.38 9.772 
19.51 20. 82 6.71 11.457 
144.70 153.71 6.22 11. 105 
41.93 44.75 6.73 25. 873 
54.68 58.91 7.74 28. 124 
193. 86 201.77 4.08 26. 559 
86. 84 94.31 8.60 43. 195 
97.74 102. 58 4.95 44.827 
243. 12 254.91 4. 85 45.311 

表 4 步 长 0.01 的 各 模 态 计算 结果 

实测 值 / ” 步 长 0.01 ”计算 误差 ”计算 时 间 / 

Hz 计算 值 /Hz % s 

15.75 16. 18 2.73 17.574 
19.51 20.20 3.54 19. 130 
144.70 149.72 3.47 19.946 
41.93 43.42 3.55 32. 483 
54.68 56.55 3.42 35. 801 
193. 86 200. 59 3.47 33.722 
86. 84 90.93 4.71 57.503 
97.74 101.47 3.82 55.390 
243. 12 249.71 2.71 57. 375 


表 5 步 长 0.005 的 各 模 态 计算 结果 


Tab.5 Calculation results of each mode with steps 0. 005 

焉 术 各 次 “实测 值 / 。 步 长 0.005 计算 误差 /计算 时 间 / 
Hz 计算 值 /Hz % s 

一 阶 挥舞 15.75 16.11 2.41 24. 643 
一 阶 摆 振 19.51 20.11 3.08 23.182 
一 阶 扭转 144.70 149.28 3.17 23.206 
二 阶 挥舞 41.93 43.35 3.39 45. 442 
二 阶 摆 振 54.68 56. 54 3.40 42. 160 
二 阶 扭转 193. 86 200. 46 3.40 46. 844 
三 阶 挥舞 86. 84 90.79 4.55 61.757 
三 阶 摆 振 97.74 100.91 3.24 69.708 
三 阶 扭转 243. 12 248.73 2.31 64.316 
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表 6 步 长 0.001 的 各 模 态 计 算 结 果 


Tab.6 Calculation results of each mode with steps 0. 001 


实测 值 /” ” 步 长 0.001 “计算 误差 / ”计算 时 间 / 


模 态 名 称 计算 信 /H 人 

一 阶 挥舞 15.75 16. 12 2.35 37. 149 
一 阶 摆 振 19.51 20. 09 2.97 39. 626 
一 阶 扭 转 144.70 149.26 3.16 37. 803 
二 阶 挥舞 41.93 43.33 3.37 67.717 
二 阶 摆 振 54. 68 56. 49 3.35 64. 465 
二 阶 扭转 193.86 200. 30 3.32 62. 683 
三 阶 挥舞 86. 84 90.76 4.53 86.291 
三 阶 摆 振 97.74 100. 88 3.18 80.090 
三 阶 扭转 ”243.12 248. 67 2.29 91.591 


表 7 步 长 0.0005 的 各 模 态 计算 结果 
Tab.7 Calculation results of each mode with steps 0.0005 


实测 值 / ” 步 长 0.0005 计算 误差 / ”计算 时 间 / 


A Hz 计算 值 /Hz % 

一 阶 挥舞 15.75 15,.12 2.35 71,332 
一 阶 摆 振 19.51 20. 09 2.97 87.644 
一 阶 扭转 144.70 149.24 3. 14 94.155 
二 阶 挥舞 41.93 43. 32 3.32 187.042 
二 阶 摆 振 54.68 56.49 3.31 173.125 
二 阶 扭转 193.86 199.75 3.04 158.004 
三 阶 挥舞 86.84 90.73 4.48 362. 711 
三 阶 摆 振 97.74 100. 85 3.18 329.415 
三 阶 扭转 243.12 248.63 2.28 291. 620 


计算 相对 误差 /% 


率 1 


条 二 时 必 格 闪 罗 


而 


| 0.01 0.001 0.000 1 
步 长 


图 3 不 同步 长 下 的 各 模 态 固有 频率 计算 误差 


Fig.3 Calculation error of natural frequency of 


each mode under different step size 

为 了 进一步 考察 计算 步 长 减 小 对 计算 精度 提高 
与 计算 时 间 代 价 的 影响 ,此 处 对 表 2 ~ 表 7 的 计算 结 
果 进 行 了 进一步 统计 ,给 出 了 步 长 减 小 对 计算 误差 
与 计算 时 间 的 影响 分 析 表 ( 见 表 8)。 表 8 中 : 
si0-1234560 代 表 6 种 不 同 的 计算 步 长 ;s, 代表 每 种 
步 长 下 前 三 阶 模 态 计算 值 与 实测 值 之 间 的 平均 计算 
误差 ;t; 代表 每 种 步 长 前 三 阶 模 态 的 平均 计算 时 间 ， 


单位 是 s;As,, = |s;,1 -s; | 代表 步 长 减少 量 ;As,， = 
1s -2i| 代 表 随 着 步 长 的 减 小 ,每 种 步 长 相对 前 
一 种 步 长 所 带 来 的 平均 计算 误差 减少 量 ; Ai,,, = 
| 二 = 天 | 代表 随 着 步 长 的 减 小 ,每 种 步 长 相对 前 一 
种 步 长 所 带 来 的 平均 计算 时 间 增 加 量 ;As As 
= |ei, -2i|/11iw -ti 代表 随 着 步 长 的 减 小 ,每 
种 步 长 相对 前 一 种 步 长 所 带 来 的 平均 计算 误差 减少 
量 与 平均 计算 时 间 增 加 量 比值 , 即 由 于 步 长 减少 所 
带 来 的 计算 精度 提高 与 所 消耗 计算 时 间 代 价 的 比 
值 ,该 比值 越 大 表明 所 采取 步 长 的 计算 收益 和 效率 
越 高 。 从 表 8 中 可 以 看 出 ,在 本 研究 选取 的 6 种 步 
长 中 ,计算 步 长 取 0.01 时 ,Ae;,1,/Ati, 最 大 , 即 步 
长 取 0.01 所 带 来 的 计算 收益 和 效率 最 高 。 
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> 8 步 长 减 小 对 计算 误差 与 计算 时 间 的 影响 分 析 


CTab.8 Analysis of the influence of step size reduction 


Cs on calculation accuracy and calculation time 
ss 
三 A girl 
| ei ti Asiri1 | 
Atirl 


0.1 7.9 9.9 

0.05 6.03 27.36 0.05 1.87 17.46 0.11 
0.01 3.49 36.56 0.04 2.54 9.2 0.28 
0.005 3:22 44.58 0.005 0.27 8.02 0.03 
0.001 3.17 63.05 0.004 0.05 18.47 0.003 
0.0005 3.12 195.01 0.0005 0.05 131.96 0.0004 


”综合 表 2 ~ 表 8 和 图 3 ~ 图 4 的 分 析 结果 ,本 研 


究 选用 0.01 作为 基本 迭代 步 长 开展 后 续 旋 站 轴 心 
运动 轨迹 计算 。 
3.2 刚性 旋 收 的 轴 心 运动 轨迹 预测 

将 旋翼 视 为 刚性 时 ,在 已 知 旋 壳 各 模 态 响应 的 
状态 下 ,下 面 将 计算 刚性 旋 吕 在 低速 状态 ( 桨 尖 速 度 
为 198. 12 m/s) 和 高 速 状态 ( 桨 人 尖 速 度 为 365.76 m/s) 
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下 的 浆 载 轴 心 轨迹 ,分 析 旋 豆浆 载 的 运动 情况 。 
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图 5 低速 状态 刚性 旋翼 (收敛 的 ) 轴 心 运 动 轨迹 
Fig.5 Axis trajectory of rigid rotor (convergent) at low speed 

图 5 中 横 纵 坐 标 分 别 表 示 效 载 轴 心 的 在 2 个 方 
向 上 相对 偏 移 量 ,计算 数据 均 已 使 用 桨 叶 弦 长 数据 
相 除 进行 量 纲 归 一 化 处 理 ( 图 6 ~ 图 8 同 )。 分 析 得 
知 , 在 低速 状态 条 件 下 ,启动 一 段 时 间 后 直升机 刚性 
旋 避 轴 心 运动 为 一 近 圆 轨迹 , 且 有 不 断 向 稳 态 点 运 
动 的 趋势 ,证 明 此 时 旋 踊 轴 心 运动 趋 于 收敛 ,机 体 结 
构 稳 定 。 


5 0 5 

* 

图 6 高速 状态 刚性 旋 经 (发散 的 ) 轴 心 运动 轨迹 

Fig.6 Axis trajectory of rigid rotor (divergent) at high speed 
如 图 6 所 示 ,在 高 速 状 态 条 件 下 ,刚性 旋 豆 轴 心 

运动 轨迹 为 一 发 散 的 近 圆 形 , 桨 载 有 偏离 起 点 的 趋 

势 , 说 明 此 时 直升机 结构 出 现 风险 ,将 产生 旋翼 的 结 

构 损 伤 ,可 能 引起 严重 事故 。 


3.3 和 柔性 旋翼 的 轴 心 运动 轨迹 预测 


将 旋 辟 视 为 柔性 时 ,在 已 知 旋翼 各 模 态 啊 应 的 
状态 下 ,下 面 将 计算 柔性 旋回 在 低速 状态 ( 效 尖 速度 
为 198.12 m/s) 和 高 速 状态 ( 桨 尖 速 度 为 365.76 m/s) 
下 的 桨 坑 轴 心 轨迹 ,分 析 旋 恤 桨 载 的 运动 情况 。 

从 图 7 可 以 看 出 ,在 低速 状态 条 件 下 ,柔性 旋 辟 
由 于 式 向 和 了 向 阻尼 刚度 不 同 , 奖 载 轴 心 的 整体 俯 
视图 呈现 出 近似 椭圆 形 , 其 半 长 轴 和 半 短 轴 随 着 时 
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间 推 移 逐 步 减 小 且 趋 于 几何 中 心 ,证 明 此 时 旋 波 轴 
心 运动 趋 于 收敛 ,机 体 结构 稳定 。 


15 


nh 
己 


10 坊 


| 


7 低速 二 状态 和 性 旋 届 ( 收 全 的 ) 轴 心 运动 轨迹 


Fig.7 Axis trajectory of flexible rotor (convergent) at low speed 
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Fi@ 3 Axis trajectory of flexible rotor (divergent) at high speed 
SN 从 图 8 可 以 看 出 ,在 高 速 状态 条 件 下 ,柔性 桨 载 
el 圆 形 , 且 极 远 点 不 
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结构 不 稳定 ,有 解体 风险 。 


1 ) 针对 Johnson 倾 转 旋 站 机 不 平衡 载 答 前 飞 动 
力学 模型 ,推导 了 直升机 旋 惨 弹性 阻尼 和 惯性 力 综 
合作 用 时 桨 叶 的 挥舞 和 摆 振 运动 方程 ,并 利用 车 加 
原理 得 到 了 效 载 轴 心 运动 方程 ,给 出 了 直升机 旋 波 
的 轴 心 运动 轨迹 求解 流程 。 

2) 以 某 型 直升机 旋 波 系统 为 例 , 验 证 了 本 研究 
所 提出 旋 中 桨 叶 模 态 分 析 方 法 的 准确 性 ;从 固有 频 
率 计算 误差 和 计算 时 间 的 角度 出 发 ,给 出 了 计算 效 
率 最 佳 的 求解 步 长 选取 方法 。 

3) 预 测 了 典型 飞行 状态 下 的 桨 载 轴 心 运动 轨迹 ， 
为 直升机 旋 渗 系统 设计 提供 基础 方法 和 评估 手段。 
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